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Abstract. 

The dynamical equations of an electromagnetic field coupled with a conducting 
material arc studied. The properties of the interaction are described by a classical 
field theory with tensorial material laws in space-time geometry. We show that 
the main features of superconducting response emerge in a natural way within the 
covariance, gauge invariance and variational formulation requirements. In particular, 
the Ginzburg-Landau theory follows straightforward from the London equations when 
fundamental symmetry properties are considered. Unconventional properties, such as 
the interaction of superconductors with electrostatic fields are naturally introduced 
in the geometric theory, at a phenomenological level. The BCS background is also 
suggested by macroscopic fingerprints of the internal symmetries. 

It is also shown that dissipative conducting behavior may be approximately treated 
in a variational framework after breaking covariance for adiabatic processes. Thus, 
nonconservative laws of interaction are formulated by a purely spatial variational 
principle, in a quasi-stationary time discretized evolution. This theory justifies a 
class of nonfunctional phenomenological principles, introduced for dealing with exotic 
conduction properties of matter [Phys. Rev. Lett. 87, 127004 (2001)]. 
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1. Introduction 

Classical electrodynamics has been a commonplace in several areas of mathematical 
physics. Thus, many of the physical problems linked to the applications of function 
theory or differential geometry are taken from electromagnetism. In particular, we recall 
the very elegant formulation of the electromagnetic (EM) field in space time geometry. 
The full set of Maxwell equations may be simply expressed as dj-" = and 5JF = J', 
for JF (the electromagnetic field tensor) a closed 2-form defined on the 4-dimensional 
Minkowski space, the current density 1-form, and 6 the so-called codifferential 
operator p. The closedness condition of (rfjF = 0) allows a local representation 
in terms of a potential 1-form A, introduced by = dA. Maxwell laws may then arise 
as the Euler-Lagrange equations for an action functional in terms of A. Within this 
geometric formulation of electromagnetism, the consideration of symmetry properties 
simplifies and illuminates the theory. In particular, covariance (which assures invariance 
of the equations under transformations of the Lorentz group) is very convenient and 
gauge invariance (invariance under transformations of the potential A ^ A + dx) must 
be satisfied. 

All the above ideas are mainly established in the study of the electromagnetic 
(EM) field in vacuum, but can also be of great help in the research of interactions 
with electric charges within macroscopic media. Thus, it is known that the most 
relevant aspects of superconductivity: expulsion of magnetic fields, zero resistance, flux 
quantization and the phase-voltage relationship at the gap between superconductors, 
may be straightforwardly obtained from gauge invariance considerations [2]. 

Recently, and motivated by some experimental puzzles in superconducting physics, 
several additional tools, conventionally restricted to the research of EM fields and sources 
in vacuum, have been applied for the description of material laws. To be specific, 
space time covariance of the phenomenological equations of superconductivity has been 
considered as the appropriate framework for explaining the still unclear interaction 
of these materials with electrostatic fields. Essentially, and inspired in the principle 
of Lorentz covariance, when the electric and magnetic fields are treated at the same 
level, one predicts both electrostatic and magnetostatic field expulsion with a common 
penetration depth A |Sl IH E]- This formulation has been used jE] as a basis for 
explaining the so-called Tao Effect, an intriguing experimental observation, in which 
superconducting microparticles aggregate into balls in the presence of electrostatic fields 
[7j. However, both theoretical objections jHj as well as experimental results P raise 
concerns on the universality of the common A treatment. 

We want to emphasize that the topic of electrostatic field expulsion was already 
addressed by the London brothers (TU] in the early days of superconductivity. 
Nonetheless, owing to the lack of experimental confirmation , they eventually decided 
to formulate their celebrated equations of superconductivity in a noncovariant form. 
Such a lack of relativistic covariance leads to theoretical difficulties, but they are 
conventionally avoided by postulating the absence of electrostatic charges and fields 
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within the sample. This point of view has been adopted by the scientific community 
over decades, until the revived interest mentioned in the previous paragraph. 

In the light of the above comments, it is apparent that finding a physically sound 
covariant expression for the material physical laws, in intrinsic geometric terms may be of 
help. In fact, it will be shown that this will not only release the theoretical dissatisfaction 
caused by using a somehow amended theory, but it will also allow to incorporate new 
physical phenomena in a natural way. In an effort to produce the most concise and 
general equations, we propose in this work, the use of differential forms on Minkowski 
space, within the area of electrical conductivity. However, it should be emphasized 
that relativistic exactitude is not meant to influence the kinetics of superconducting 
carriers. What one tries to do is to write the physics in the clearest form so as to get 
information on the underlying mechanism, through the macroscopic (phenomenological) 
electrodynamics. Being electrodynamics a fundamental interaction in the case of 
superconductors, one tries to scrutinize within the nature of the phenomenon, taking 
advantage of the EM field symmetries. 

As a first approximation to the problem, we will restrict ourselves to linear effects 
in the material response, i.e. nonelectromagnetic degrees of freedom enter through 
the linear laws = Q ■ T oi J = S ■ A + ut, with O G T^^ and S G field- 
independent tensors, and u a gauge adjusting 1-form. Although simple, these ansatzs 
provide celebrated material laws. 

The study of the law J' = E ■ A + u, with the requirement of gauge independence, 
will straightforwardly lead to the phenomenological London ^2] and Ginzburg-Landau 
(GL) equations in covariant form, for the trivial case SJ;^ = aS'j^ (here, 5^ stands for the 
Kronecker's symbol). This will serve as a basis for postulating a Lagrangian density, 
still a controversial topic related to the time dependent GL theory ^2]- On the other 
hand, nonequal diagonal terms in H will suggest relativistic BCS effects, in the manner 
introduced in [Hj. Unconventional degrees of freedom, such as electrostatic charges will 
be unveiled in our gauge independent proposal. 

The study of the law J = Q ■ T will lead to the concept of dissipative interaction, 
and thus to the absence of a direct variational formulation. However, we show that 
under quasistationary evolution, and following prescribed covariance breaking, one may 
issue a 3D variational principle, in terms of differential forms, at least for certain 
nonconservative interactions. This concept will also allow us to address some open 
questions in superconductivity. In particular, we justify the use of restricted variational 
principles in applied superconductivity for the so-called hard materials. At a first 
approximation, these superconductors are treated by a nonfunctional law, in which 
the electrical resistivity jumps from zero to infinity, when a critical value in the current 
density is reached p!^ . 

The work is organized as follows. First, some mathematical background material, 
regarding notation and operations with differential forms is recalled in section |21 The 
presentation is conceived so as to provide the minimal tools for using the powerful 
geometric language in what follows. Also, the variational formulation of pure or coupled 
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electromagnetic problems is reviewed. For further application, we emphasize the relation 
between gauge invariance and the structure of the Lagrangian density. Section El is 
devoted to the proposal of covariant and gauge invariant material laws for conducting 
media [J^iJ-") and J' (A)]. We show that a variational formulation for the covariant 
Ohm's law does not exist within the electromagnetic sector, and that, on the contrary, 
superconducting dynamics finds a very appropriate host in such formalism. Then, 
in section |3] we give the rules for breaking covariance in the formulation from the 
geometrical point of view. We introduce a restricted variational theory and exploit 
the benefit of using it by applying these ideas to exotic superconducting systems, in 
which the material law is a nonfunctional relation. The implications of our analysis and 
further proposals for the theoretical studies on conducting materials are summarized in 
section El 

Rationalized Lorentz-Heaviside units will be used for the electromagnetic quantities 
along this article (with c = 1). 

2. Mathematical background 

2.1. Covariant formulation of electromagnetism 

Covariance is a well established requirement for fundamental physical theories. In the 
more general statement it means that the physical laws can be formulated in intrinsic 
geometric terms, and therefore they are independent of local coordinate descriptions, 
the reference system. There are sometimes further requirements associated to some Lie 
group invariance, e.g., the Galilean group in Classical Mechanics or the Lorentz group 
in special relativity. 

2.1.1. Basic notation and definitions The classical geometric treatment of the Maxwell 
equations is formulated in a space-time manifold M. endowed with a fiat Lorentzian 
metric g [here we choose the signature (— , +, +, +)]. Recall that a pseudo-Riemannian 
metric g m a, manifold M provides us with an isomorphism of the C°°(M)-module of 
vector fields in that of 1-forms, g : X(M) /\^{M), given by {g{X),Y) = giX,Y), 
which allows us to transport the scalar product from X{M) to /\^(M). We shall use the 
same notation g for such a product 

5((ai,a2) = 5'(?~^("i),?"^("2)) , ai, 02 G /\ ^(M). (1) 
We also recall that given an n-dimensional pseudo-Riemannian manifold M, we define 
for any index p < n a. C°°(M)-linear map star from the C°°(M)-module of p- forms into 
that of {n — p)-forms by means of (see e.g. |14j ) 

aAi.(3={-iyg{a,l3)Vo\ Va,/3G/\P(M) (2) 

where g{a, (3) is the scalar product in f\ ^(M), i.e. if a and (3 are decomposable p-forms, 
a = ai A ■ ■ ■ A Op and /3 = /3i A ■ ■ ■ A /3p, then g{a, j3) = det{g{ai, (3j)). Vol is our choice 
of the volume form associated to the metric; it is given by 

Vol = ± V(-l)" det ^ rfx^ A ■ ■ ■ A dx"" , (3) 
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where s denotes the signature, the number of negative squares appearing in the quadratic 
form associated with g when written in its diagonal form. 

In affine coordinates for the particular case of Minkowskian space A4 (for which 
n = 4 and s = 1), {x'^ | /i = 0,1,2,3}, the signature is 1 and for a = a^dx^ and 
P = Pfj. dx^ in /\^{A4), we find (recall summation over repeated indices) 

g{a, (3) = g^^a^Pi, = -aoPo + aiPi (4) 

and we make the choice 



Vol = ^/-detg dx^ A dx^ A dx^ A dx^ = rfx^ A dx"^ A dx^ A dx° , (5) 

so that in a local basis 

7^ (rfx°) = -dx^ A dx"^ A dx^ , * (dx^) = -dx^ A dx^ A dx^ , 
-k{dx^) = dx^ A dx^ Adx^ , * (rfx^) = A dx^ A dx^ , (6) 

and more generally 

i,(3 = - p^dx^ A dx^ A dx"^ - (3i dx^ A dx^ A c?x° + /32 dx"^ A A 

- /^s A c/x^ A dx"^ . (7) 
Similarly for 2- and 3-forms, the star operator is determined by 

7^ (rfx° A dx^) = -dx^ A dx^ , ★(rfx° A dx^) = dx^ A dx^ , 

*(cix° A dx^) = -dx^ A cix^ , i^{dx^ A (ix^) = dx° A dx^ , 

^(rfx^ A dx^) = -dx^ A rfx^ , -k{dx^ A c/x^) = dx° A rfx^ , (8) 

and 

* (rfx^ A rfx^ A dx^) = -c?x° , ^ ((ix° A c/x^ A dx"^) = -dx^ , 

*(dx° A dx^ A dx^) = c^x^ * (rfx° A rfx^ A dx^) = -dx^ . (9) 

By defining the pairing for scalars as (?(!, 1) = 1 we also find 

^(1) = rfx^ Arfx^ Arfx^ Arfx°, *(dx° A dx^ A dx^ A dx^) = -1 . (10) 

Notice that if G /\ ^(M), then {-k o 'k)uj = (— 1)*'("~p)+'* a; and, in particular, in the 
Minkowskian case for which n = 4, s = 1, {-k o 'k)uj = (— l)^+^Ci;. 

As usual, the exterior derivative of a p-form over an ra- dimensional manifold is 
defined by 

du = ?^^!iLdL dx^ A dx'' A ■ ■ ■ A dx'" . (11) 
dx^ ^ ' 

On the other hand, the codifferential operator is defined by 5 = -kd-k^ it maps p-forms 

into {j) — l)-forms, and do b ^ b o d'vs, the Laplacian (or D'Alembertian) of the manifold. 

The space of the physical variables T for the EM Classical Field Theory is the 

set of closed 2-forms in M.. Let us introduce some geometric objects and manifolds 

appropriate for a detailed description of the EM theory. This will serve as a link between 

the geometric language and the more conventional (analytical coordinate dependent) 

statement of the problem. It is well known that 2-forms in M. are sections for the 
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vector bundle A : T*Ai A T*A4 ^ A4 of skew-symmetric (0,2) tensors over M. For 
a given fibre bundle 'j : E ^ M, we denote by J^7 the first jet bundle J3] of 7, the 
manifold of equivalence classes of sections s : M ^ E with first degree contact at a fixed 
point m G M. J'^'-f is the natural geometric framework for a system of first order partial 
derivative equations (PDE). Elements of J^7 are denoted by j^s, the equivalence class 
of all sections with zero and first order partial derivatives at the point m equal to those of 
s. For a given section s, j^s denotes de first jet lift, a section j^s : M — > J^7 defined by 
j^s{m) = jl^s. Local coordinates in E adapted to the projection 7, {x'^,y"'} with {x'^} 
coordinates in M, determine associated local coordinates in J^7, {x'^,y°',z'^ = d^t^y"'}- 
Now, it is more transparent that a system of first order PDE H{x'^,y"',dxt^y"') = 
represents just a submanifold P C J^7. 

In particular, local adapted coordinates for A : T*M. A T*A4 — A4, {x^, J-'^u} with 
jF(m) = T^u{m)dx^ A dx" , determine local coordinates for J^A 

Closed 2-forms, i.e., sections T for A fulfilling dT = 0, are such that their first 
jet bundle lift j^jF : Ai J^X goes to zero through skew-symmetrization. The 
skew-symmetrization in J^A is a natural map sk-^ : J^A (T*Ai)^^ defined by 
sk^{il^J-') = dT{m). In local coordinates, 

sk^ o j^T{m) = i (^^,,, + g,^,, + g,^,^) {m)dx'' A dx"" A dx" . (13) 

Similarly, the codifferential S determines a natural map r : J^A ^ T*Ai given by 
'^Um^) = 5J-'{m). A local coordinate expression can be obtained through the previously 
presented set of coordinate relations for the star operator. 

2.1.2. Maxwell equations Below, we show that the standard expressions of the Maxwell 
equations in terms of vector calculus operators are recovered from the previous formalism 
when a particular coordinate system is specified. The set of Maxwell equations is a 
system of first order PDE in the EM field JF. Then, they can be geometrically described 
as a proper submanifold P C J^A, or equivalently, as a family of geometric tensorial 
equations in J^A whose set of solution points determines P. More precisely. Maxwell 
equations are just dJF = and 5J^ — = 0, where JT" represents a the 4-current 
density, either prescribed or related to JF through some material law. If J' is prescribed 
J : M ^ T*M, then P can be alternatively defined as P = sA;3"^(0) fl r~^(lmj'). 
On the other hand, if there is some relation J{T), J : T*M A T*M T*M (not 
necessarily linear), then P = {j^jF|s/c3(j^JF) = and T{jl^J^) = j7'(JF(m.))}. 
In local affine coordinates 

= -Ei dx° A dx' + ^eijkB' dx^ A dx'' (14) 

determines the electric and magnetic vector components of JF, obviously coordinate (i.e. 
reference frame) dependent. Here eijk is the totally skew-symmetric, Levi-Civita, tensor. 
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Both = and 5jF = J' represent the covariant version of Maxwell equations. 
Let us obtain their coordinate dependent version. 
The equation dj-" = is given by 

dT = - d^jEi dx^ A dx' A dx^ + ^e^j-fc d^oB"- dx° A dx^ A dx'' 

+ ^(^ijk d^'B\ dx^ A (ix^ Adx'- = 0. (15) 
Using the vector differential operator V = {d^i, d^^, d^a), dj-" = becomes 

(16) 



V- B = 

V X E + d^oB = 0. 



On the other hand, 5JF = ^ is given by 

5J^= - {d^iEi + 9^.2^2 + d^aE^)dx^ + {-d,^sB2 + d^2B^ - d^oEi)dx^ 
+ (-9,i53 + d,,Bi - d^oE2)dx^ + (-9,25i + 4i52 - d^oE^)dx'' 
= J = -pdx^ + Jidx' (17) 

and, in vector analysis notation, 

VxB-9,oE = J 

V.E = p. (^^^ 

Notice that the chosen metric tensor (— , +, +, +) determines the 1-form representation 
for J', so that its corresponding 4- vector field is ^{J^) = pd^o + g^^Jjd^i. On the other 
hand, recall that 6J^ = 0, the continuity equation {dtp + V ■ 3 = in standard notation), 
is a consistency requirement for 5jF = JT", easily obtained from 6^a = {—lY'^^-kd'^'ka = 
for any p-form a G /\^. For further use, we give below the wave propagation expressions 
both in intrinsic terms and in standard notation. They arise by taking the exterior 
derivative in equation (fTTj) . i.e. 

ddJ^ = dJ. (19) 

This admits the following coordinate form in terms of the D'Alembertian operator for 
our metric (□ = —d^^ + V^) 

°^ = -^^J (20) 

Let us now review the geometrical notation and properties of the vector potential. 
The closedness condition djF = can be (locally) integrated, and disappears from the 
theory, by describing the EM field as the exterior differential of a 1-form A, the 4- 
potential 1-form, JF = dA. The potentials are local sections for the cotangent bundle 
IT : T*Ai Ai, and JF = dA is determined by the composition of the first jet lift 
and skew-symmetrization map, JF = sk2 o j^A, with sk2 '■ J^i^ T*M. A T*Ai. We 
have natural adapted coordinates {x'^, A^^, A^^^^, = dx-^A^^} at J^tt, and {x'^,Af^,F^u = 
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A^ i, — A^^^} at T*Ai A T*M.. Then, the electric and magnetic vectors are represented 
through the potential by 



with A = —0 dx^ + Ai dx\ Note that in terms of our chosen metric tensor (— , +, +, +) 
the vector field corresponding to the potential 1-form is ^{A) = (pd^o + g^^ Aj d^^. 

On the other hand, the local representation of the EM field through a local integral 
potential 1-form gives way to nonuniqueness, with local gauge invariance symmetry 
A' = A + dx] quotient by this gauge invariance determines the classical physical degrees 
of freedom. Thus, it is important to notice that the potential A is not directly a physical 
quantity at the classical level. However, recall that topological obstructions (JF being 
closed but not exact) associated to the configuration of the problem have sometimes 
experimental consequences; this is, for instance, the case of the well-know Aharonov- 
Bohm effect, in which charged quantum particles are influenced by the circulation of A in 
a multiply-connected region where the electromagnetic fleld vanishes. From the physical 
point of view, gauge invariance is an additional requirement of the theory whenever 
the potential appears, either in the material laws (phenomenological or fundamental), 
interaction of the EM fleld with conducting samples, or in the variational theories, where 
the action functional must be gauge invariant. 

2.2. Variational principles of electromagnetism 

The existence of a variational formulation is well known for fundamental theories 
and mandatory if one wishes to connect the classical and quantum levels. In many 
cases, additional degrees of freedom interacting with the EM fleld will also have 
their own kinetic and potential Lagrangian terms, and the energy conservation law 
reflects the possible transfer between the EM and other energies. However, it must be 
stressed that phenomenological theories do not always permit a variational formulation, 
because dissipation of the EM energy can be balanced by generation of other kind of 
energy (usually thermodynamical) , and the corresponding degrees of freedom may be 
disregarded in the theory, i.e., the system under consideration is open. This section 
is devoted to review some examples of Lagrangian theories for the EM fleld: free, 
interacting with prescribed 4-current, and with an additional scalar fleld. A number 
of speciflc features will be outlined for their application to the proposal of material laws 
in the following one. 

2.2.1. Free EM theory The Lagrangian function for the free Field Theory, L : J^tt M, 
is given by 



Notice that, being in fact a real function in T*M. A T*A1, L o .sk2 is constant along 
the flbres of sk2 '■ .7^7 — > T*M. A T*Ai. This property is nothing but the gauge 



E = - 
B = V X A 




(21) 



L Vol = A = ^ {E^ - B^) Vol 



(22) 
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invariance at an algebraic level, L takes the same value for two jets j^Ai and j^A2 whose 
difference j^A2 —j^^Ai (the jet bundle has a natural affine structure) is symmetric, i.e., 
sk2{jlnAi) = sk2{jlnA2)- This also means that the Lagrangian is singular, and given a 
solution section A for the Euler-Lagrange equation, we can build physically equivalent 
solutions by adding to A arbitrary sections whose first lift is in the kernel of sk2- 
Obviously, are nothing but closed one forms. Maxwell equations are first order PDE 
in the EM field JF, and therefore are second order PDE in the potential A. Euler- 
Lagrange equations for a first order Lagrangian are second order. The second order jet 
bundle J^7 of a fibre bundle {E, 7, M) is defined in a similar way to J^7, taking now 
equivalence classes of sections up to second order derivatives. Now, second order PDE 
are in geometric terms submanifolds of J^7. The geometric description of the Euler- 
Lagrange equations for a classical field theory goes through the Poincare-Cartan form 
0£ = dsylj + LVol (see [12] for a detailed description of the geometric treatment of 
the Euler-Lagrange equations in classical field theories), which is an n-form in J^7 for 
n = dimM, and Sy is a vector valued n-form generalizing the vertical endomorphism in 
tangent bundles, defined for each volume form V in M. From O^, the Euler-Lagrange 
form, S-^ = (72)*((iL A Vol) + dhQ]^ with dh the total derivative mapping r-forms in J^7 
into (r + l)-forms in J^7, happens to be an (n + l)-form in J^7. The geometric equation 
{j'^s)*[Si2 = for unknown s represents the second order partial differential equations 
(PDE) fulfilled by sections s of 7 : i? — > M making stationary the action functional 
S{s) = /L(j^s)V, that is, they are the Euler-Lagrange equations of the variational 
principle. In local coordinates {x^ , y"" , z"^, z'^,^} in J^7 

(f)} rl f^} \ d 

9y« dxi'dz^J ^ dxi" " ^ - ^ ^ 

and its components are the well known Euler-Lagrange equations in coordinate form. 
For the case of EM field theory, we find Gauss' and Ampere's law in vacuum 



(24) 



and 



+ S''eukf-'""''r%, (9,™ A„) = (25) 

On the other hand, a direct application of Noether's theorem, connected with the 
invariance of the action under the Lorentz group, leads to the concept of canonical 
energy-momentum tensor 
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As usual, a symmetrized and gauge invariant version is preferred, in order to ease 
interpretation. Thus, we will use the field symmetrizing technique fH] 

T^" = Q^g'^'' + , (27) 

from which the conservation law dT^'^/dx^ = follows immediately. 

The 00 component of the (symmetrized) energy-momentum tensor T^^" ^H] is 
= i(E2 + B2), the classical EM energy. In local coordinates, the zero component of 
dx'^T'^'^ represents the conservation of EM energy, and the spatial components are the 
conservation of EM momentum. In particular, the continuity equation 

(9^oT°° + V ■ (E X B) = (28) 

is the balance between the energy density time variation and power flow; integration of 
E X B along the boundary of a compact region measures the power transfer out of it. 
Obviously, this is not a conserved quantity for non Lorentz-invariant Lagrangians. 

2.2.2. Electromagnetic field with external sources In the presence of sources (prescribed 
electrical charge and current densities) for the EM field, an additional term A A -kj' in 
the Lagrangian, i.e. 

LVol = A - A A i^J (29) 

determines the modified Euler-Lagrange equations 

6T = J. (30) 

Notice that LVol = A i^J^ ~ A A i<J is no longer constant along the fibres of the 
skew-symmetric projection sk2. Thus, for a gauge transformation A^ A + dx, 

LVol ^ LVol + dxAi<J = LVol + d{xi^ J) - xd* J ■ (31) 

The current density must fulfill the continuity equation 6J^ = in order to maintain 
the gauge invariance of the action functional. Then, the Lagrangian is modified by a 
divergence term, which does not affect the dynamics of the system. We stress that gauge 
invariance is a fundamental ingredient of the theory, not only forcing consistency for the 
Maxwell equations, but also determining transformation properties for the Lagrangian 
densities or the material laws. 

In local coordinates the term —A A -kj' takes the form 

(-p0 + J-A)Vol (32) 

as it can be easily computed from the action of * presented above. 
In this case, the Euler-Lagrange equations take the form 



and 



(9T 

+ S'^euke^'^'^dx. (a,..AO = _ = , (34) 
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obviously equivalent to equation (fTH|). 

Eventually, the energy balance equation becomes 

9^oT°° + V ■ (E X B) = -E ■ J , (35) 

showing that there is a transfer between EM energy and other modes. This 
may correspond to reversible storage of energy by charged particles, irreversible 
thermodynamical losses, etc. E ■ J is sometimes called thermodynamical activity. 

2.2.3. Coupling with a scalar field: the Klein-Gordon equation In some problems, the 
current density will not be prescribed, but arise as a consequence of charged particles 
moving in the EM field according to the electromagnetic force. New degrees of freedom 
have to be incorporated to the Lagrangian density by writing J in terms of the particle 
positions and velocities, and adding a kinetic term for the masses of the particles and 
possibly a potential interaction term between them. If, at a macroscopic level, the 
number of particles allows to consider a continuum charge and current density, the 
total system will be described by the EM and fluid fields. This could be a good 
approximation for some physical systems such as low density plasmas. Nevertheless, 
the usual interaction of EM with matter is still unsatisfactorily described in this way. 
Charged particles (electrons) move in a material lattice, with which they interact, and 
interchange momentum and energy. Thus, new degrees of freedom (vibrational, for 
instance) should be incorporated to the model. This may cause serious difficulties, and 
in some instances a phenomenological material law may be of great help. 

Just as a starting point for our subsequent proposal (section l3.3j) . we recall the 
simplest theory that couples the electromagnetic phenomenon and a relativistic material 
field ip in covariant form. Thus, if one considers a scalar spin-0 field representation for 
the dynamics of the charged particles, the basic Lagrangian ^7] may be written as 

£Vol= -J^A^J^+^di>Ai.dip-—ijijVol-V{i;i;)Vo\ + Cintyo\ (36) 
2 Sttiq 2 



with 



■ - - g2 _ ■ 

iq(ip dip — ip dip) — —ipipA 

h 



(37) 



2mo 

Note that C has been split up as Cem + + indicating the self-interactions of 
both fields, and a coupling term between them. The potential term V{il)ip) may also 
incorporate reversible interactions of the particle field with other degrees of freedom, as 
the underlying material lattice. 

As one can easily check, the Euler-Lagrange equations of the corresponding action 
integral, when A, ip and its complex conjugate ip are taken as independent variables, are 
nothing but the coupled Maxwell and Klein-Gordon equations. Just, one has to identify 
the electromagnetic current density (take dj^ in the Lagrangian) with the KG current 
density times the basic electric charge (g), plus a vector potential related term, i.e. 



2mn 



2q' - 

iqiip dip — Ip dip) —ipipA 

n 



(38) 
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3. Covariant material laws in conducting media 

Interaction of the EM field with matter is described by some material law, an equation 
determining the response of the medium through the appearance of charge and current 
densities under an applied EM field, usually generated by sources which can be typically 
considered far away form the area of interest. In principle, in order to propose simple 
material laws one must take care of preserving the basic rules of any EM theory, that 
is, covariance and gauge invariance. As already discussed, the existence of a variational 
principle depends on the possibility to treat additional degrees of freedom associated to 
other kinds of energy that balances the possible dissipation (or creation) of EM energy. 
In general, this will not be the case for a phenomenological theory in which we are 
exclusively interested in the EM sector, but there should be a Lagrangian density for 
more fundamental theories that try to consider all the degrees of freedom present in the 
system under study. 

Below, we present a material law fulfilling explicitly the first two requirements, 
covariance and gauge invariance, which could describe the classical response of 
conducting matter, but which does not allow a variational principle. Afterwords, we 
will introduce several laws allowing a variational formulation, which will immediately 
lead to the concept of superconductivity. The geometrical treatment will allow a natural 
upgrading of the theory, so as to infer the covariant GL equations, as well as a first 
indication of the BCS background. 

3.1. Nonvariational tensorial laws J[T^ 

In the following, we will consider an external EM field Te generated by some given sources 
outside a particular region Q of the space, which are incorporated to the problem through 
appropriate boundary conditions on the boundary dQ^ and a sample material in i?o C Q, 
generating an additional EM field Tr as a response to the applied excitation. Here dQ 
is taken far away from i?0; so that the EM material response J> can be neglected there. 
The local current density J within the sample, with compact support in i?o; will be 
determined by some material law J[T\ T = Tr- Therefore, we consider that the 
physical system is not isolated, being fed by the external sources through the boundary, 
and possibly with dissipation in the sample, i.e. transfer of EM energy into another kind 
(thermodynamical, mechanical, etc.). In general, the Lagrangian density obtained by 
replacing J in the Lagrangian used for prescribed currents will not produce the correct 
Euler-Lagrangian equations. In fact, widely used material laws determine Maxwell 
equations which are not variational. Then, dissipative force densities are added to the 
free Euler-Lagrange equations by just writing the current density as J{T^ in equation 
(^H))- The simplest choice is a local first order approximation. In order to preserve the 
geometric flavor, a quite general law (under the pointed restrictions) may be handled 
as the tensorial equation, by introducing a (2, 1) type tensor 6 contracted with the EM 
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(0, 2) tensor JF to determine the 1-form current density JT": 

txSj^={J',txe) = txJ ^ J^ = e;^j-,p (39) 

Different components of the G tensor represent well known EM versus matter interaction 
behaviors. Thus 

-p = Ql'J^OJ + &'^J'^, (40) 
contains both electric and magnetic polarizability, while 

J^ = e^^ok + ef (41) 

describes Ohm's law and magneto-conductivity. 

In passing, we note that the celebrated covariant form of Ohm's law ^H] J'^ + 
u^u^Jij = {1/ R)J^^'^Uy, where is the fluid four- velocity, corresponds to the previous 
equation when = 0, m° = 1 and O^"* = (1/-R) 5/. 

Eventually, the set of (Maxwell) equations to be solved in this obviously covariant 
and gauge invariant model of electromagnetic interaction are 

c/J^ = 0, ix5^= (^,^xe) VX. (42) 

For instance, if one goes again to the case of linear isotropic conducting media, in the 
absence of electrostatic charges, the wave diffusion equations d5T = d{Q ■ JF) become 

V^E = \ 

Rdt dt^ 

(43) 

V^B = \ . 

R dt dt^ 

They are the well known equations describing the penetration of electromagnetic fields 
in conducting media, and have to be solved supplemented by boundary conditions for 
the fields. 



3.2. Variational scalar laws J^{A) 

Let us now consider a less classical material law, determined by a covariant relation 
between the current density and the local potential field, J^{A). It is apparent that an 
additional current u will have to be considered in order to preserve gauge invariance of J^. 
In principle, such additional current may look a purely mathematical artifact; however, 
as we will see, widely accepted classical models of superconductivity are obtained in 
this way, and the new current can be understood as the classical shadow of a more 
fundamental quantum dynamical superconductivity theory. This geometric approach 
is therefore a natural way to generate classical approximate models for macroscopic 
quantum properties of EM interaction with matter. 

We will consider purely variational problems, so that EM energy variations are 
balanced by the energy variation of the additional field. At a first stage, we will analyze 
some simple choices of the additional current to be incorporated to the interaction term, 
and will not consider the origin of this new current, i.e., the underlying field and its 
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corresponding kinetic and self-interaction terms will be neglected in the Lagrangian 
dynamics. Later on, the new field will be incorporated to the theory by a minimal 
coupling prescription. 

3.2.1. The London model (dui = 0) Let us consider first the following EM Lagrangian 
density, with a very simple choice of interaction term, 

L = A - A , (44) 
Z Za 

where is a 1-form 

J = aA + uj, (45) 

with a a constant parameter, and uo a closed 1-form [du = 0) associated to some field 
not discussed yet. The EM Euler-Lagrange equations become 

6J^={aA + uj)=J, (46) 

determining Maxwell equations for the proposed material law. In local coordinates, with 

A = —(f) dx^ + Ai dx^ and u; = cuq dx^ + ujidx\ we have 

d.^,Ei = {acj) - uq) , 
(V X B)i = d^oEi + {aA, + 

In order to preserve gauge invariance for J', when one considers the gauge 
transformation A^ A' = A + dx, should transform according to 

uj uj' = uj — adx ■ (48) 
Under this rule, the Lagrangian is manifestly gauge invariant. 

Topological properties The integral of (spatial part) along a curve inside the sample, 
has two components, coming from the potential aA and the new current u. Notice that 
a A + u; is equivalent to a{A + dx) + {uj — adx)- Then, the circulation of JT" contains 
two gauge invariant components. In the present model, the circulation of uj vanishes for 
trivial topologies because of Stokes' theorem. 

Note that, the new current being closed, there is a gauge fixing where uj locally 
vanishes, by an adequate selection of x- However, the u;— independent theory, which may 
be identified with the basic London equations jTHj is somehow unsatisfactory, because 
one loses physical information. In particular, notice that, although locally vanishing, uj 
can however contain a topological charge for non trivial topologies, like a hole on a plane 
as configuration space (or an infinite cylinder in 3-D space), i.e. 

rfcj = 0^w~V/^ (l uj-d\ = 0. (49) 
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Continuity At this stage, the continuity equation, SJ' = 0, is not a consequence of 
gauge invariance, because we have not considered yet the full action functional. Here 
6J^ = is a consistency condition for Maxwell equations, and determines the relation 
a6A + 6uj = between the EM potential and the new current. According to this, if one 
takes the cj— independent formulation, the gauge fixing freedom is lost, and one should 
work within the Lorenz gauge condition 5 A = 0. 



Wave equations By using the definition of JF = dA and applying the exterior derivative 
to the material Maxwell equations, in order to eliminate the new current, we get 

d5J^ = aT . (50) 

The left hand side represents de D'Alembertian of the EM field [6dT = 0), so that 
we have obtained a wave propagation equation with sources inside the sample. In 
local coordinates, and after appropriate identifications of the parameter a and the flux 
expulsion length scale A one has 

B d^B 



A2 dt^ ' 

(51) 

E d'^E 



A2 df" 

Recall that the gauge-independent wave equations are insensitive to the 1— form uj, that 
disappears by the closedness condition. As E and B are observable quantities, these 
equations are a test for the soundness of the model in which u is closed 

Under quasi- stationary experimental conditions, where the wave propagation can 
be disregarded, we get the celebrated London equation 

V^B = g. (52) 

Additionally 

V^E=^ (53) 

represents a penetration of electric field, not usually considered, but mandatory for 
covariant considerations [31 IH El ■ 

Although u has been eliminated in the current model, it becomes clear that it is 
an unavoidable geometric ingredient to deduce London's equations from a variational 
principle while maintaining the gauge invariance. In the following section, we study the 
more general case in which du 0. A different physical scenario will arise. 



3.2.2. The modified London model (du ^ 0) A more general model is obtained by 
relaxing the closedness condition for u, i.e., we allow for du ^ 0. Now, one can find 
non vanishing circulations for u even for trivial topologies. Solutions for the material 
Maxwell equations 

6J^ = aA + uj = J , (54) 
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plus the continuity equation 

S{aA + uj) = (55) 
will determine particular distributions of usual EM and the full current density. 

Topology The contribution of u to the circulation relates to magnetic flux inside the 
sample (recall that ^ A ■ d\ = JjB ■ ds). This represents the existence of vortices 
(quantized magnetic flux) associated to the so-called type II superconductors [12]. The 
spatial distribution of u into compact regions where du ^ and a surrounding space 
with duj = becomes a simplified classical model for the existence of quantum vortices 
in type II superconductors. The fact that circulations embracing the vortex regions do 
not vanish is of importance, as it carries information about the superconducting field 
related to uj. 



Wave equation The wave equation reads d6J^ = aT + duj^ i.e. 



B S^B 



V X a; , 



at 



(56) 



Recall that the static or quasi-static approximations are obtained by neglecting either all 
the time derivatives or just the second order ones, in the previous formulas for covariant 
superconductivity. It is important to notice that in any case, both V x a; and Vcjq 
should be maintained, generalizing previous proposals jH]. Also, recall that though uj 
has been introduced for mathematical consistency, its temporal and spatial components 
[uj^^iS) become observable charge and current densities. By the moment, they have 
to be considered as phenomenological quantities, but below they will acquire a (more 
fundamental) microscopic significance. 



3.3. A covariant, gauge invariant and variational model of superconductivity [J'{A,tp)] 

Let us now look for a field such that its associated current density u fulfills the former 
requirements of material law for EM interaction with matter. As we have already 
shown in section El the complex Klein-Gordon field ijj represents a simple choice for 
that purpose, because it fulfills covariance, gauge invariance and variational formulation 
requirements. The above introduced free parameter a will be adjusted so as to identify 
the correct gauge transformation rule, and the interpretation of ipip as the density 
of superconducting carriers will bring us to the famous Ginzburg-Landau equation, 
which here is proposed in a covariant and gauge invariant framework. Recall that the 
superconducting carriers (Cooper pairs) are spin-0 combinations of electrons, and thus, 
the KG equation seems to be a reasonable starting point for a covariant field theory 
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of superconductivity. Thus, identifying ip as a. charged KG field, the London current 
density term u may be expected to be 

uj = - — dip — ip dip) . (57) 
Outstandingly, upon gauge transformations, u verifies the required law 

UJ ^ uj' = UJ + '^dx (58) 
if one defines the rules 

^ ^ Q-^ix/f^^ ■ ^ ^ e'l^/^ip (59) 

and 

1 - ^^^^ (a(^\ 
^ = = -a . (60) 

y mo 

To this point, mo and q are just an effective mass and charge for the KG particles. The 
transformation rules in equation ()59|) correspond to the internal U{1) symmetry of the 
charged field. 

As it has been discussed elsewhere [21 E] , the KG particles may be interpreted as 
mediating Higgs bosons, whose fingerprint in the theory is the mass term A A -kA for 
the electromagnetic potential. In conclusion, the generalization of the covariant London 
Lagrangian should read 

CVo\= ljr/\^jr+ J_^^/\^/)^_V^(^^)Vol-^7/;7/;Vol, (61) 

as it follows from a minimal coupling principle, applied to the field ip. To be specific, 
we have used the covariant derivative 

D = d + i^A (62) 

Now, taking variations in equation ()6H1 respective to the 1-form A, and to the 
scalar fields ip and ip produces the set of Euler-Lagrange equations 

d dC _ dC 



d dC _ d£ 

dxf^ d{dip/dxt^) dip 
t 

d dC _ dC 

dxf^ d{dip/dxf^) dip 
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h 



2mo 
t 



+ i^A^ 
n 



d^-i^A^ 
n 



2mo - 

When the self interaction model is chosen to be 



- dV mo - 



(63) 



(64) 



one gets the covariant Ginzburg-Landau-Higgs equations of superconductivity In 
geometric notation we get 



(65) 



for the Maxwell equations {J' includes a term due to the redistribution of carriers 
density), and 



for the carriers wave function in superconducting state. 



(66) 



3.3.1. Energy-momentum Starting from equation (j6T| and keeping in mind that C 
depends on the fields A'^,ip,il!, one may calculate the full energy- momentum tensor 
from the expression 



dC OA, 



dC dijj dC dtjj 



+ 



(67) 



which will account for the self-energy of the electromagnetic and KG fields, plus the 
interaction. The calculation results in the symmetrized form 



2^-2 



2mo 



q - 
doipdoijj + -^Aoipi/j 



(6^ 



This is nothing but the GL free energy with relativistic effects. Recall that the 
nonrelativistic conventional expressions have been augmented not only by the particles 
rest energy, but also by the electrostatic terms E"^ and 

2 



Aoij 



A V/2 



(69) 



as one could obtain from equation pijl in the London limit. 



3.3.2. Nonrelativistic limit A straightforward calculation allows to obtain the 
nonrelativistic limit of the GL Lagrangian. This will produce a Schrodinger-like equation 
for obtaining the low frequency limit of the time-dependent Ginzburg-Landau (TDGL) 
theory. 
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As a starting point, we split up the wave function in the form 

7/;(x, t) = $(x, t)e-^™o'='*/^ (70) 
with $ the nonrelativistic part of the wave function, for which the relations 

<moC^|$| , |gAo$| < moC^|$| (71) 

must hold. Notice that in the previous formulas c is not normalized, so as to ease 
quantitative comparison. Equations (j71|) mean that, compared to the rest energy, the 
nonrelativistic energy may be neglected, and that the potential is flat enough so as to 
avoid spontaneous pair creation. 

By starting with equation (pT|l . separating the /i = and fi = 1,2,3 components, 
and implementing the above relations, one obtains 

Cnonrel ^ -\f^uF^^ - ^ ($9o$ " ^8,^) + -^D<l ■ D<1> - V($<l) , (72) 

4 2 ^ ' IniQ 

with D the spatial part of the covariant derivative. We remark that this Lagrangian 
produces the TDGL equations for the nonrelativistic limit, as well as the correct limit for 
time-independent solutions, in which case one recovers the conventional GL free energy. 
This point may be easily checked, just by examining the Euler-Lagrange equations. 
When ones considers variations respect to $, the familiar Schrodinger-like equation for 
$ is obtained. 

9$ h"^ / ia 
ih— = - — V + ^A <1> + ^<1> + I.|<1>|'<1> (73) 
at ZrriQ \ fi / 

Recall that the inclusion of the mass term for the relativistic KG field has been essential 
in order to reach the above result. 



3.4- Tensorial material law J{p ■ A) 

The wave equations (jSH), or their quasi- stationary approximations and 
determine the penetration profiles for the electromagnetic field in a Type I 
superconductor. As mentioned before, equation is the widely accepted London 
model for the penetration of the magnetic field and its electric counterpart, which 
has been disregarded for decades since the pioneering works of the London brothers 
[Tin IH] ! who concluded that it wasn't physically sound to consider a finite decay length 
for the electric field. 

Following the geometric orientation of this paper, we state that maintaining the 
magnetic penetration depth and rejecting the electric one (making it zero, as an ansatz) 
is inconsistent with the idea of covariance. This has also been recalled in references 
[SI in El Ej , where several theoretical proposals along the lines of our study are given, 
some of them also connected with the experimental reality [3 [H]. However, a new 
controversy has appeared in the literature about the subject of an electric penetration 
depth in superconductors. There are firmly grounded theoretical reasons, as well as new 
experiments (see jH]) which support that, though nonzero, the typical penetration length 
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of the electric field is negligible with regard to the typical one for the magnetic field. 
Roughly speaking, non superconducting charges in the material, which do not play a 
role in the equilibrium electric current transport, must be taken into account when an 
electric field generates an electrostatic response in the sample. The order of magnitude 
of the E— penetration depth related to the normal charge density distribution happens 
to be much smaller than its magnetostatic counterpart, and this explains the lack of 
evidence for the electrostatic phenomenon, at least within the experimental conditions 
considered. 

In our phenomenological approach to the subject, the former physical considerations 
can be incorporated to the material law for the superconducting state in a covariant and 
gauge invariant way. As a preliminary proposal, let us consider an a bit more general 
material law than the one introduced in section IT2l To be specific, let us concentrate 
on a still linear, but tensorial relationship J^{A) of the form 

J = E-A + Q (74) 

with the (1,1) tensor S depending upon two phenomenological constants in the form 
= a, Sq = a + P and vanishing non diagonal components, in a suitable affine 
coordinate system (the rest frame for the sample). The intrinsic tensor S could be 
computed in arbitrary coordinate systems through the standard transformation law for 
(1, 1) tensors. 

Gauge invariance Notice that an additional 1-form Q has been added to the material 
law for gauge invariance considerations. In particular, under a gauge transformation 
^ 1-^ ^ + dx, the corresponding gauge transformation for Q becomes 

n^n-E-dx- (75) 

Within the above model, we have the Maxwell equations 

dj^ = 6J^ = J = E- A + n. (76) 

As stated before, these should be completed with the continuity equation 6J' = 
5(S ■ A) + 6Q = 0, which gives information about the properties of the additional 
current Q. On the other hand, the simplifying hypothesis du = assumed in the first 
London model (section 1^.2.11) cannot be translated into this model because 

dn = and d{n-E- dx) = (77) 

will be, in general, inconsistent for an arbitrary gauge function i-e., one can have 
d{E ■ dx) ^ for a nontrivial tensor S. 



Wave equations A straightforward computation allows to obtain the particular form 
of the wave equations (jlOI) for this case. One has 



□ B = — -VxO (7^ 
A"' 
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and 

with the definitions a = — and l3 = —\jv^. 

We emphasize that, under a static configuration (all the time derivatives are 
neglected), the former equations become 

V^B = — - V X (80) 

and 

V^E= (l + l^E-Vfio^^E-Vf^o. (81) 



As proposed in jH] , the electrostatic penetration depth combines the effects of the 
magnetic one (A) and of another phenomenological constant [y). The microscopic origin 
of such dichotomy has been discussed in that work within the BCS theory framework. 
At the level of this article, what we can state is that is the macroscopic manifestation 
of internal degrees of freedom beyond the Ginzburg-Landau f/(l) gauge invariant model. 
Just note that the tensorial gauge transformation rule in equation ()75|) does not allow 
to introduce a complex scalar field ensuring a gauge invariant theory in the manner of 
equations ()58j) and . Further aspects of this problem will be discussed elsewhere. 



4. Noncovariant material laws in conducting media 

In the previous section, we have exploited the concept of relativistic covariance for 
studying electromagnetic material laws in conducting media. Having settled the basis 
for the geometrical description in 4-dimensional Minkowski space, we are ready to 
discuss about appropriate restrictions to 3-dimensional Euclidean space. This will be 
done below. The motivation for this part is to justify the use of restricted variational 
principles in the study of quasistationary conduction problems. 



^.i. Breaking space-time covariance: reference frame 

In this section, a particular reference (rest or laboratory) frame will be chosen, allowing 
to decompose the EM 2-form JF into its electric and magnetic vector field components. 
From a geometric point of view, a reference frame corresponds to the choice of local 
coordinates (preferably affine) in the space-time manifold Ai. However, in order to 
maintain the freedom about the spatial coordinates (spatial covariance), we will consider 
the splitting Ai = 'Rx Q, where Q is the three dimensional Euclidean space, or an open 
submanifold with boundary if the particular properties of the system make it desirable. 
M represents the absolute time for the laboratory frame, and we have both natural 
projections tti and 7r2 of M x Q over each factor. Also, by fixing a time value t, there 
is a trivial morphism jt : Q ^ A4, jtiq) = (t, q), allowing to pull-back forms in A4 into 
forms in Q. By doing it for each point of an interval [0, T], we can define one-parameter 
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families of forms in Q. Free selection of spatial frames means that the EM theory is 
now developed in the 3-D tensor analysis framework. When necessary, we will also fix 
a gauge in order to simplify some equations, but gauge invariance of the theory must 
be maintained even after breaking covariance. Different components of the space-time 
forms are obtained by puUing-back the original form (getting the space-like component) 
or its contraction with T = dt (time-like component), the natural vector field on the 112 
fibres M. 

With this notation we can obtain the corresponding components of the EM field, 
the potential and the current density, associated to the reference frame. The Hodge 
operator in Q, with the Euclidean metric qe = Sij (Kronecker's delta), will be denoted 
by * in order to distinguish it from the space-time * Hodge operator in Ai, and the 
exterior differential in Q will be denoted by d instead of the d in Ai. Similarly, the 
(Q, 5'£)-codifferential will be denoted by S. 

Now, we can define the vectorial and scalar fields in Q 

j7(^) = *B(t), j7(zr^) = -E(t), 
j:{A) = A{t), j;(zr^) = -0(t), 

j;(J)=J(t), j7(zrjr) = -p(t), (82) 

with the obvious identification of magnetic and electric vector fields, vector and scalar 
potential fields, as well as electric vector current and charge densities. Spatial r- 
dependence of the fields has been avoided in the previous definitions for simplicity. 
Maxwell equations in the reference frame take the geometric form 

dF = — > {(5B = 0, dE + 9f(*B) = 0} , 

6J^ = J — y {SE = p, *dB-9tE = J}. (83) 

Quasistatic limits, with some field constant in time can be considered in the previous 
equations. 

If one neglects c?t(*B), i.e., electromagnetic energy is only stored in electric form, 
one reaches the so-called EQS (ElectroQuasiStatic) regime, in which 

(5B = , dE = 

SE = p, *dB -dtE = J. (84) 

On the contrary, the MQS (MagnetoQuasiStatic) approximation corresponds to 
neglecting dtE. Then 

(5B = 0, dE + 9j(*B) = 

SE = p, *dB = J. (85) 

Apparently, the EQS and MQS regimes arise when some characteristic speed in the 
problem is small as compared to c. In the general case, when no speed is neglected, 
energy is alternatively stored either in electric or magnetic forms, and one has a 
propagating wave. In the next section, we will concentrate on systems for which the 
MQS limit is attained. Our scenario will be as follows: some initial magnetostatic 
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configuration is perturbed, giving place to a transient process in which electric fields 
and possible charge densities appear. Then, the system is driven to a final magnetostatic 
configuration, and remains there until perturbed again. Dissipation J ■ E can appear in 
the transient process and, although small, it cannot be neglected in the study. 

4-2. Spatial variational principles in quasistationary processes: the law J{T^ 

As it was seen before (sections 12.21 and Ej), the genuine variational formulation of 
electrodynamics is done in and in terms of the potential 1-form A. However, some 
physical systems are successfully analyzed in fixed reference (laboratory) frames, while 
keeping spatial covariance. For instance, this is trivially true when one focuses on the 
static equilibrium configuration of conservative systems. Then, minimization of energy 
produces equations determining the fields E and B. What we show below is that such 
idea may be generalized to the quasi- stationary evolution of dissipative systems. Under 
certain conditions, dynamical equations produced by spatial variational principles are 
justified. 

In Classical Mechanics it is well known that an initially conservative system which 
is slowly drifted by an additional small non conservative force, linear in the velocity, 
admits an approximate variational principle. This represents an adiabatic evolution, 
in which the energy, although not conserved, varies slowly according to the adiabatic 
parameter. More specifically, let us consider the dynamical equation 

with small parameter A. By adding to the Lagrangian L = {l/2m){dx/dtY — V{x) the 
so-called Rayleigh dissipation function j^HI {l/2)tX{dx/dt)'^ we find the Euler-Lagrange 
equation 

d^x „ ^ ^ ^ dx , d'^x , 

which differs from the correct one in a negligible term for small time intervals, because 
both t and A are considered small. 

Within a time interval [0, AT] the increments fulfill the equation 

mA(^) + {d,V)AT = -XAx - ]-ATA{A . (88) 
dt 2 dt 

The adiabatic hypothesis can be reformulated by saying that the dissipative force 
Xdx/dt varies slowly along the evolution. In Classical Mechanics, this is often used for 
conservative systems with periodic orbits in which an adiabatic evolution generates small 
variations of the parameters on each cycle, e.g., the Poincare map. It can also be used 
to perform a numerical integration through time discretization and minimization of the 
approximated action functional at each step, allowing to apply minimizing techniques for 
the integral, usually more reliable that a direct numerical integration of the differential 
equations. Recall that the Euler-Lagrange equations are just stationarity conditions for 
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the action integral, but in any reasonable physical variational system, the solutions are 
local minimizers for the action, although possibly not global ones. 

The above method can also be performed for an EM system with slow drift 
between stationary configurations through a transient material response to small source 
variations. The procedure could be developed for quite general systems but, in order to 
fix the ideas and taking into account the particular application which follows, we will 
consider an MQS approximation for a (type II super-)conductor, with vanishing E and 
p in the initial and final configurations. The transient evolution will generate a small 
electric field, with possible local charge density production, that will be neglected. 

We denote by Bq and Jo the initial stationary values, obviously fulfilling V x Bq = 
Jq. Bo(r) and Jo(i') fields are known functions, and we do not need to introduce a 
potential field for the stationary configuration. A potential vector field A is chosen 
so as to describe the transient evolution, while the scalar field is ignored in this 
approximation and should be determined by the null charge density prescription. This 
can be interpreted as the selection of the temporal (Weyl) gauge. Then, along the 
evolution one has 

B(t,r) =Bo(r) + AtB(r) with AtB(r) = AB(t, r) 

J(t,r) = Jo(r) + AtJ(r) with AtJ(r) = AJ(t, r) . (89) 

Additionally 

AtB = V X A, E = -dtA (90) 

as the representation in terms of the potential. 

Faraday's law is automatically verified (geometric equation) while Ampere's law is 
the one to be (approximately) determined through a variational principle. The system 
under consideration will fulfill the following conditions: 

1) the finite sample material occupies the region Rq, Rq d Q with Q also finite but 
dQ far away form Rq so that the material EM response decays to zero in dQ, and the 
monitored sources are out of Q determining the feeding of the system through boundary 
conditions of the magnetic field in dQ. Such experimental conditions are those of a PDE 
control type problem, with a vectorial distributed parameter B^ form the sources and 
control dynamical equations. Maxwell equations, determining the response of the system. 

2) the adiabatic hypothesis determines slow variations of the sources form an initial 
value Bso to a final value Bsi in a time interval [0, T], with corresponding slow evolution 
Bq ^ Bi and Jq ^ Ji. Correspondingly, the transient electric field E(t) is also small. 
The final values are determined by the variables to be used in the analysis A^B and At J, 
with Bi = Bq + Aj^B, Ji = Jq + A^J. We also neglect the EM wave by considering 
instantaneous material response, that is, with typical time response negligible with 
regard to the typical time parameter of the control variable. Moreover, and similarly to 
the previous mechanical example, along the adiabatic evolution the electric field, playing 
the role of the dissipative term, can be considered constant, i.e., E = E(r) in [0,T] and 
V X B = J. 
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The Maxwell equations must be complemented with some material law, that we 
consider in the form A J = G(E). With E being small along the time interval, we 
make the hypothesis that we can properly approximate the material law by linearizing 
it to AJ ^ K - E, with K a. 3 x 3 matrix representing the Jacobian of G in the origin. 
Inhomogeneous K{r) could be considered, but here we will choose the case of K constant 
within the sample and vanishing outside. 

Inspired by the above mentioned mechanical Lagrangian, we define an MQS 
Lagrangian density for (M x T*Q,7t,'R x Q) in order to determine a variational field 
theory for the system 

L = X + hdtA'' ■ K ■ dtA (91) 

where the super-index T denotes the transpose. The associated Euler-Lagrange 
equations become 

V X (V X A) + 9t itK.dtA) = , (92) 
and working on them by substitutions we get 

V X AtB = K -E + tdtiK -E) ^ AtJ . (93) 

Above dfE is neglected by the adiabatic hypothesis. Notice that we have obtained 
Ampere's law, while Faraday's law was already fulfilled by the potential representation. 

The next step is to perform the time integration along [0, T] in order to get a purely 
spatial principle. We have chosen T small for a better approximation when neglecting 
tdt{K ■ E). This allows to perform an approximate integral by just considering mean 
values according to the initial and final values of the magnetic field and current density, 
as well as constant E within the interval. Writing the Lagrangian in terms of these 
variables we find the minimization principle 

mmjvolj^ (^^{AtBY + ^tE.AtJ^dt (94) 

After integrating in time according to the above prescriptions, and avoiding global 
numerical factors, we get 

min / ((AtB)2 + T < E > -A^J) vol (95) 
Jq 

a purely spatial principle. 

For the sake of completeness and consistency, let us check that, under the hypothesis 
considered, the spatial variational principle reproduces the correct dynamics. For further 
application, we will do that in an unconventional way, rewriting the spatial Lagrangian 
in terms of the variable A^B by using both Ampere's law and the linearized material 
law. Notice that, contrary to the prescription of Faraday's law through the use of the 
potential in the space-time variational problem, we are here prescribing Ampere's law 
by using it in the substitutions of the spatial Lagrangian. 
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Starting from 




(96) 



and performing the Euler-Lagrange equations for the field theory in the spatial variables, 
we obtain 



that is, the discretized version of Faraday's law. 

We emphasize that the above property is mainly grounded in the time discretization, 
where the vector potential A can be rewritten as — TE, an integration which cannot be 
directly fulfilled in space-time, and in the fact that the scalar potential has disappeared 
form the formulation. In fact, the reader can write the spatial Lagrangian in terms of 
the vector potential A, i.e., prescribing the time discretized Faraday's law, and obtain 
in a more conventional way Ampere's law as the Euler-Lagrange equations for this 
Lagrangian. We have sketched this in ()91fl^nj) 

For a more physical interpretation of the result in equation ()95|) . one can identify 
a magnetic inertial term and a dissipative term, which are balanced in order minimize 
the addition of magnetic flux changes and entropy production [T^ . 

4-3. Application to hard superconductivity: variational statement for nonfunctional 
{E,J} laws 

Some physical systems are better described by general relations (graphs) between 
their variables, rather than by functional ones. As a particular case of technological 
interest, we recall the conduction property of the so-called hard type-II superconductors. 
According to the phenomenological Bean's model j2j|, the scalar components of E and 
J, when currents flow along a definite direction, are related by the nonfunctional relation 
depicted in figure H Notice that, if the electric field is nonzero at some point, one has 

J = -/cp (99) 

at such point. However, if E = any value J G [—Jc, Jc] is allowed. From the physical 
point of view, the superconductor reacts with a maximal current density flow to the 
application of electric fields. When the excitation is canceled {E — * 0) the flow may 
remain as a persistent nondissipative current. Notice that the hard superconducting 
material displays the conventional quasistatic zero resistivity, until a certain level of 
current transport is demanded (Jc). Current densities above this threshold are no longer 
carried by supercharges, and a high electrical resistance is observed. 

A fundamental justification of the above model (so-called critical state model), 
the physical interpretation of the material parameter J^. (critical current), and more 
sophisticated versions may be found in fSl 1221 12^] and references therein. 




(97) 



which, through (V x A^^B) = A^J and K ^ ■ A^J = E, becomes 




(98) 
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Figure 1. {E, J} graph (conduction law) for a hard type-II superconductor, according 
to Bean's model. Vertical lines correspond to infinite resistivity above Jc- 



Here, we recall that variational methods are especially useful for solving and 
generalizing the above statement. To start with, we interpret it as a more realistic 
limiting process (see figure^). Thus, the hard material allows lossless subcritical current 
flow, while it reacts with a high resistivity E = R{J — Jc) for J > (analogously for 
negative values of J). The harder the material, the higher value of R, and the more 
realistic the graph law approximation. Then, one can start with equation (jHSI) and 
notice that, as R becomes larger, the second term also increases with J > J^. In the 
limit of infinite slope, this fact can be taken into account by reformulating the variational 
principle as 



da 



-Jr 



da 



+Jc 




Figure 2. Restriction sets for the current density (J G ri), corresponding to the 
behavior of hard type-II superconductors. Optimal control solutions obey the condition 
J* e dn. 
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mill S = 




for |V X B| < Jc . (100) 

The inequality (unilateral constraint) determines that the mathematical framework for 
the model is the so-called Optimal Control Theory [25 , an extension of the classical 
variational calculus for bounded parameter regions. In the Optimal Control language, 
we have a performance (cost) functional S to be minimized under the control equation 
V X B = J for bounded parameter |J| < Jc- As a very relevant property of this 
variational interpretation, we remark that the control region for the parameter may be 
understood as a physically meaningful concept. Thus, one may pose the problem in the 
very general form 



min S = 




for V X B G 1] , (101) 

with r2 C M'^ some restriction set, prescribed by the underlying physical mechanisms. 
The conventional statement, given by the graph in figure Q is nothing but the particular 
case Q = [— Jc, Jc]- This is depicted in figure |21 In the literature, several possibilities 
for the set Q have been studied, and identified as the fingerprint of different physical 
properties. For instance, elliptic restriction sets have been shown to reproduce 
experimental observations related to anisotropic current flow [22]; a rectangular set 
has been used for producing the so-called double critical state model j^Hl, in which two 
independent critical current parameters (parallel and normal) are used, etc. 

Technical note The use of the principle (IHKjl as a basis for obtaining (jlOlj) deserves 
some explanation. Not all the hypotheses used in the former case are straightforwardly 
translated. In particular, one has to recall that the linearization AJ = K ■ E has to 
be considered for excursions of J around J^. However, according to figure ^ when the 
electric field is reversed one has A J ~ 2 Jc. Again, the finite jump may be smoothed by 
taking a small time step. Then, the size of the region where A J = 2 Jc is negligible, and 
the fault has a very small weight in the integral to be minimized. 

Finally, we recall that control equations linear in the parameter, as the current case 
of interest, produce the so-called bang-bang solutions characterized by the condition 

J* e on (102) 

i.e. the optimal solutions take values at the boundary of the allowed set. For instance, 
the control variable jumps between 1 and —1 when Q = [—1, 1] (see figure |21). 

Example: the infinite slab in parallel field. The optimal control approach to critical 
state problems in superconductors has been sometimes misunderstood and qualified as 
a more or less intuitive approach of restricted applicability. Among other questions. 
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it has been said that it lacks information, because the electric field is absent of the 
theory. However, this quantity is in the essence of the variational statement, which has 
been obtained from the general Lagrangian, under the quasi-stationary assumptions 
f section l4.2j) . Below, we show with an example that, indeed, the variational statement 
contains the electric field. Thus, when the minimization process is performed, one 
obtains the condition that J belongs to the boundary, and also a set of Lagrange 
multipliers (associated momenta in the Hamiltonian formalism) closely related to the 
electric field. The comparison of our equations with a more conventional approach that 
directly uses an E(J) relation will prove this aspect. 

Let us consider the infinite superconducting slab depicted in figure El < 
d/2] \y\, \z\ < oo), and the problem of determining the electromagnetic response to 
an excitation field parallel to the surface. Owing to the symmetry, one may use the 
assumptions 

B = (0,5,(x),i?,(x)) 

J = (0,J,(x),J,(x)). ^ ^ 



Assuming isotropic conditions, equation (jlOlll takes the form 

min S =- I (AB)\ol 
2 Jo 



for J 



OB 



eQ=[0, Jc] . (104) 



dx 

Recall that the cost function depends on the field increment AB = B — Bq for each step 
of time, in the evolution of the system. 

According to Pontryagin's maximum principle [21], this problem is solved by 
combination of (i) the canonical equations for the associated Hamiltonian 

H = p-AJ-i(AB)2, (105) 

and (ii) AJ* = J* - Jq such that H{A3*) > H{AJ) , \/ J E n, i.e., 

C?B D 

max[p- AJ] ^ — = Jc-- (106) 

ox p 



[B^(x), B^(x)] 



z© 



'fx 



Figure 3. Superconducting slab, subject to magnetic field parallel to the surface. 
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This leads to the system 

dBy ^jPy 

dx ^ p 
dB^ _ jPz 
dx p 

dx 

^^AB. (107) 

By using the definition {py,Pz) = ^t{Ez, —Ey) the system may be rewritten as 
dBy _ J Ez 

dx E 
dB^^ _jEy 

dx E 
dE,, AR, 



dx At 

(108) 

dx At ^ ^ 

Taking derivatives, and inserting the standard notation of dots and primes, one obtains 

92 R E,E - E,E 

= Jr— — 

dxdt E^ 

d^ By ^ E,,E — E,iE 



and 



dxdt E'^ 



j^ ^y^ (109) 



and thus. 



dxdt 
pa r> 

EyE EyE 



E^ 

4'=^^^^^^. (Ill) 

Now, using polar coordinates in the plain, 

Ey = E cos (p 

E, = Esmip (112) 

and taking space and time derivatives, it follows 

E" = E" cos Lp — 2E'ip' sin (p — E((p'Y cos ip — E^p" sin 
Ey = E" sin ip + 2E'(p' cos (p - E^f sin ip> + E^" cos if 
Ey = E COS (p — E(p sin (p 

Ez = E simp + E(p cos (f . (113) 
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Eventually, back-substitution in equation (jllip leads to the system 
E" = [^'?E 

These differential equations, together with suitable boundary conditions, allow us 
to obtain the penetration profiles for the modulus of the electric field and for its angular 
direction. They have been obtained from our variational approach, and fully coincide 
with the expressions obtained by the E(J) method in reference |27j . 



5. Conclusions and outlook 



In this article, we have shown that the geometric formulation framework of the classical 
electromagnetic field in terms of differential forms in Minkowski space may be extended 
to the study of conducting materials. 

Related to recent intriguing experimental observations, and to the underlying 
dissatisfaction caused by a noncovariant theory for several aspects of superconducting 
electrodynamics, we have proposed a new phenomenological approach to the problem. 
We show that with the geometric jargon, the theory may be highly simplified. In the 
language of 1-forms, covariant superconductivity is merely a linear law which admits 
the inclusion of phenomenological constants and physical quantities. Such quantities 
allow a direct physical interpretation, as they are a part of wave equations in which 
they couple to the observable macroscopic fields E and B. 

Having clarified the basis of a covariant theory, we also present the complementary 
side of how covariance should be broken if required by the mathematical counterpart of 
some physical process. Thus, we show that quasistationary conduction problems may 
be treated in a spatial 3D covariant framework by pullback of the Minkowski space 
1-forms to the Euclidean space. Taking advantage of this prescription, we have been 
able to justify the use of restricted variational principles in some problems of interest 
for applied superconductivity. 

Two different possibilities for the material law have been considered, J[T^ and 
J{^A)^ i.e. the current density 1-form either depends on the electromagnetic field 2- 
form or on the potential 1-form. Being the simplest choice, linear dependencies have 
been considered. 

The linear law J{T) is obviously covariant and gauge invariant by construction. 
However, it is not variational. Only after breaking covariance, and under a quasi- 
stationary approximation, one can issue a restricted variational principle for such a 
case. In particular, we obtain an approximated spatial variational statement for linearly 
dissipating systems (nonconservative forces are proportional to 9jA, whose temporal 
variation may be considered small). 

The linear law J {A) is covariant and allows a variational statement by endowing the 
electromagnetic field Lagrangian with an interaction term of the form J /\-kJ . However, 
in this case, gauge invariance has to be required. By doing it, one is naturally led to add 
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new currents, generating the equations of superconductivity {J' = aA + u). The most 
relevant features of this phenomenon are direct consequences of internal symmetries 
in the field theory. Within the so-called London approximation, the main field is 
A. However, the 1-form u is required by the theory, coming from an additional field. 
Outstandingly, this field has observable consequences (as the presence of electrostatic 
charges and the flux quantization condition), and is sensible to the topology of the 
material. Just a step further produces the so-called covariant Ginzburg-Landau theory. 
If one identifies c^; as a Klein-Gordon like probability current density {u (x ip dip — ip dip), 
one has a conduction theory with two fields {A,ip), which may be readily identified as 
the covariant and gauge invariant generalization of the (non covariant) GL equations of 
superconductivity. Here ip represents the wave function of the superconducting carriers. 

In the first step, providing the simplest possible covariant expression for the 
conductivity of a material, we have proposed = aA + u. This fits many experimental 
observations, including the influence of electrostatic fields on superconductors [7j. 
However, accounting for other classical [T^ and very recent experiments [H], such 
proposal has been generalized to J' = S ■ A + Q, with H a (1, 1) tensor. This form 
allows to unify the referred manifestations of superconductivity by either equal or 
nonequal phenomenological constants in the diagonal terms of S. In this sense, when H is 
nontrivial, we argue that internal symmetries of the charge carriers, and gauge invariance 
are only compatible through a BCS approach. In this case, an additional field must be 
introduced, as not only the superconducting carriers are relevant. Nonsuperconducting 
charges may contribute to the static response of the material and their associated fields 
could be a matter of further research. 

Finally, we stress that the variational interpretation of a priori nonvariational 
conducting material laws under adiabatic approximation has provided us with a method 
to treat exotic materials, in which the relation between the fields is well determined 
through a graph. In particular, this has noticeable consequences in the phenomenological 
theory of type-II superconductors. We have shown that the so-called Bean's model 
for hard type-II superconductors admits a variational formulation, grounded in basic 
properties of the electromagnetic Lagrangian. This property is of utter importance in the 
field of applied superconductivity as it allows us to introduce numerical implementations 
for realistic systems, affected by finite size effects. At the level developed in this 
work, the material properties are just included by augmenting the basic term JF A 
with a dissipation function contribution. Extensions of the theory in which the base 
Lagrangian includes the conservative terms of superconductivity {J^ A-kJ', Di" ■ D$, . . .) 
are expeditious. 

Variational methods are shown to be equivalent to alternative treatments of the 
problem, but offer a number of advantages. New mathematical tools, as the optimal 
control theory, useful for discussing about the existence and form of the electromagnetic 
problem solution, as well as for hosting numerical implementations, are incorporated. 
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